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I. INTRODUCTION 

The study of supersymmetry has its origins in the 
field of particle physics, where it allows the construc- 
tion of quantum field theories containing both bosonic 
and fermionic excitations which are related by the super- 
symmetry. Such field theories are considered promising 
candidates for extensions of the standard model of par- 
ticle physics and may solve some of the open problems 
encountered in the field today. 

At the same time, supersymmetry has found appli- 
cations in many other fields of physics. In this paper, 
we will study a quantum mechanical lattice model for 
interacting fermions, which is supersymmetric by con- 
struction. The microscopic supersymmetry of the model 
makes it amenable to a variety of analytic approaches. 
Most importantly, a lower bound for the number of 
ground states can be obtained with the Witten index, 
and the exact number can be obtained from cohomology. 

The model we study was first introduced in Refs. 1 and 
2, where the model on the chain was solved explicitly us- 
ing the Bethe ansatz. Exploiting the additional tools 
that are available due to supersymmetry, further aspects 
of these models were understood: Using the Witten index 
and cohomology arguments, the number of ground states 
and some of their properties on several two-dimensional 
lattices were obtained in Refs. 3 and 4. This was used to 
demonstrate that the models exhibit surprising proper- 
ties such as superfrustration, an exponential degeneracy 
of the ground state which is not lifted by quantum fluc- 
tuations. The case of the chain was analyzed in more 
detail in Ref. 5, where the authors explore properties of 
the ground state for finite systems. The spectrum and 
its relation to superconformal field theory was explored 
in detail in Ref. 6. Recent advances include the study of 



perturbations of the model which preserve supersymme- 
try, namely staggered interactions. Interesting exten- 
sions of the supersymmetric model revealing relations to 
various spin chains were also explored. 2 ' 11 

In this paper, we will study the model on the square 
ladder. This model turns out to pose a particularly in- 
teresting numerical challenge. It has been shown to be 
critical using the behavior under boundary twists. 1 It has 
been conjectured that its continuum theory is the second 
Af = 2 superconformal minimal model, 1 but a reliable 
numerical confirmation of this conjecture has not been 
achieved so far due to strong finite-size effects in exact 
diagonalization calculations. In general, finite-size effects 
are a result of the presence of irrelevant and marginally 
irrelevant operators in the UV theory. Marginally irrele- 
vant operators are particularly notorious as they lead to 
corrections that are suppressed only logarithmically in 
the system size. One might expect that fine-tuning the 
microscopic model to be supersymmetric would prevent 
such problems from arising. We will see, however, that 
this is not the case, and that the supersymmetry imposed 
on the lattice still allows for the presence of a marginal 
operator. We resolve this problem by realising that finite- 
size corrections are strongly suppressed in the entangle- 
ment properties. We study the entanglement properties 
of the system at the supersymmetric point and a num- 
ber of perturbations away from this point, which allow 
us to establish the phase diagram and identify the super- 
symmetric point as the multicritical point in that phase 
diagram. We thus confirm numerically that the contin- 
uum limit of the supersymmetric model on the square 
ladder is described by the second superconformal mini- 
mal model. 

The outline of this paper is as follows: In Section II, 
we review the construction of the supersymmetric lattice 
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model for hard-core fermions with a focus on the rele- 
vant properties of the model on the square ladder. In 
Section III, we review the conjectured continuum the- 
ory for this model and briefly discuss the origin of the 
large finite-size corrections. This is discussed in more 
detail in Section VIII. In Section IV, we discuss the two- 
dimensional phase diagram surrounding the supersym- 
metric model, where it appears as a particular multicrit- 
ical point. 

In Section V, we turn our attention to the numerical 
methods that will be used to identify the continuum the- 
ory and study the phase diagram. Our numerical results 
for the supersymmetric point and the surrounding phase 
diagram are shown in Section VI; results for a partic- 
ular perturbation to another supersymmetric model are 
shown in Set. VII. 

In Section VIII, we discuss in some more detail the 
marginal operators present in our model, followed by a 
microscopic motivation of the field theory in Section IX. 
In Appendix A, we discuss previous work on exactly sol- 
uble limits of the model. 



II. THE MODEL 

We begin by reviewing the construction of the super- 
symmetric model of Rcfs. and and discuss some im- 
portant properties of the model. 

To construct an j\f = 2 supersymmetric model, we pro- 
vide two nilpotent operators, referred to as supercharges 
and denoted as Q and Q'. In our example, we con- 
struct these from fermionic operators that create "hard- 
core fermions" , i.e. fermions with an infinite nearest- 
neighbor repulsion. These can be written in terms of 
normal fermionic creation and annihilation operators on 
the lattice cj and Cj using a projection operator 
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where the j run over nearest neighbors of site i. Note 
that Pi commutes with Cj and cj if either i = j or i and 
j are not nearest neighbors; furthermore, Pf = Pi. The 
hard-core fermion operators are then given by 

di=PiCi 4= p i c l ( 2 ) 

We define the supercharges 

& = E C * P * Q = Y, c *P- (3) 

i i 

It is easily checked that these are nilpotent, Q 2 = 
(Qt)2 = o, and hence {Q, Q} = {Q\Q^} = 0. 

The supersymmetric Hamiltonian is now constructed 
as H = {Q\Q}- It follows immediately that [H, Q] = 
[H. Q^] = 0, i.e. the supercharges are conserved quanti- 
ties. 



A short calculation shows that 

H = E (McjP, + p^hp) + E p- ( 4 ) 

The first term is a hopping term dressed with the hard- 
core projection to ensure that only states in the allowed 
subspace are generated. The second term can be inter- 
preted as a potential term whose precise structure de- 
pends on the lattice. The Hamiltonian preserves the 
number of fermions, [H, JV rii] = 0, where rii = ctcj. 
We identify states with an odd number of particles as 
"fermionic" states, and "bosonic" otherwise. 

Supersymmetry strongly constrains the spectrum of 
the model. It is easily shown that all energy eigenval- 
ues E n obey E n > 0. In addition, the eigenvectors cor- 
responding to non-zero eigenvalues can be grouped in 
pairs , Q\ip)) , where the two states differ in occupa- 
tion by one fermion. This is reminiscent of the symmetry 
between bosonic and fermionic excitations in supersym- 
metric field theories. 

It is convenient to write the model on the square 
ladder in terms of operators d\ j (rfjj), which create 
a fermion on the upper (lower) site of the z'th rung, 
and the corresponding number operators j = d\ j d it j 

and rii j = dt.ii^j. In terms of these operators, the 
non-trivial interaction terms of the potential part of the 
Hamiltonian can be written as: 

H v2 = 2^ (n iiS n i+1 ,j +^5^4.^) (5) 

i 

H v2 > = K,l n i+2,1 + n i,in i+2A ) (6) 

i 

H v3 = ~E ( n *.J n i+l,! n ^+2,J + n i,! n l +U n l +2,l) ( 7 ) 
i 

where H v2 is a next-nearest neighbor repulsion across the 
diagonals of the plaquettes, H v2 > is a next-nearest neigh- 
bor repulsion on the same chain, and H v3 is an attractive 
three-body term. 

Many insights into the behavior of the system can be 
gained by studying the behavior under perturbations. It 
is often favorable to study perturbations that preserve 
the supersymmetry since the same analytical tools can 
be applied. A large class of such perturbations is the 
staggering of the supercharges. The latter can be modi- 
fied to have a site-dependent complex factor, i.e. 

Q f = E^ c ^ Q = E^ P ' ( 8 ) 

i i 

where the yi are complex numbers, and y~i indicates com- 
plex conjugation. None of the algebraic properties of the 
model are changed under this modification. In particular, 
the number of ground states remains unchanged for arbi- 
trary staggering. The behavior of other one-dimensional 
realizations of the model under staggering was investi- 
gated in Refs. 8 and 9. 
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The ground state degeneracy depends on the bound- 
ary conditions and the length of the system. For open 
and periodic boundary conditions the number of zero en- 
ergy states can be determined analytically. ■ For peri- 
odic boundary conditions there are 3 zero energy states 
if the length is a multiple of 4, and 1 otherwise. For 
open boundary conditions there is a unique zero energy 
state both for even and odd length of the ladder. For 
antiperiodic boundary conditions, the ground state is 
unique when the length is a multiple of 4. In all cases 
the ground state is at quarter filling. 
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Yes 


Vo,±i 


1/2 


No 




1 


No 




3/2 


Yes 




1/4 


No 


PLO-rV 0} ± 3 /2 


5/4 


Yes 



TABLE I. Relevant operators of the continuum theory. The 
second column lists their respective dimension, and the third 
column indicates whether they preserve supersymmetry. 



III. THE SECOND Af = 2 SUPERSYMMETRIC 
MINIMAL MODEL 

It was first conjectured in Ref. that the continuum 
theory that describes the supersymmctric model on the 
two-leg ladder is the second model in the series of su- 
perconformal minimal models with Af = 2 supersymme- 
try. Further support for this conjecture was presented 
in Refs. 4, 12, and 13. Let us briefly comment on the 
continuum model, before we give a short review of these 
ideas and result. 

Conformal field theories are characterized most impor- 
tantly by their central charge c and the scaling dimen- 
sions h. Without supersymmetry, the unitary rational 
CFTs with c < 1 form a discrete set referred to as the 
minimal models. It can be enumerated by an integer 
to > 3; the central charge and scaling dimensions are 
given by 



c = 1 - 
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to(to + 1) 
r(m + 1) — sto] 2 — 1 
4to(to + 1) 



(9a) 
(9b) 



where r, s are integers with 1 < r < to— 1 and 1 < s < to. 
For example, the two-dimensional classical Ising model at 
criticality corresponds to the first minimal model to = 3. 

With Af = 2 supersymmetry, a similar list of minimal 
models can be given for central charges 1 < c < 3. 15 The 
minimal series in the case of Af = 2 supersymmetry is 
also enumerated by an integer k > 1 and has 



c = 3 



6 



fc + 2 



(10a) 



p{p + 2)-r{r-2)-Ara + 2k{\-a) 2 

h = tt, ^ (10b) 

4(fc + 2) v ' 

for < p < k and r = —p, —p + 2, . . . ,p, and a = 
(a = 1/2) for the Ramond (Neveu-Schwarz) sector. In 
the Ramond and Neveu-Schwarz sectors the fermionic 
fields obey anti-periodic and periodic boundary condi- 
tions on the plane, respectively. The converse is true 
for the cylinder. The Ramond and Neveu-Schwarz sector 
correspond to periodic and anti-periodic boundary con- 
ditions in the lattice model, respectively. A characteristic 



feature of supersymmetric theories is the Witten index, 
which is related to the number of Ramond vacua or zero 
energy states. The Witten index for the k-th minimal 
model is W = k + 1. 

The theory that has been conjectured to describe the 
continuum theory of the supersymmetric model on the 
square ladder is the k = 2 supersymmetric minimal 
model with c = 3/2. An overview of other possible 
SCFTs with c = 3/2 is given in Ref. 16. Physically, 
this minimal model can be understood as the product of 
a compactified boson and an Ising theory. The operators 
in the Ising sector are 



spin (TL,h a = 1/16 
energy ip L , = 1/2 



(11a) 
(lib) 



for the left-movers and similarly for the right-movers. 
From the left- and right-moving bosonic fields, &l,r> we 
form vertex operators 



V„. 



c ,i(m+n)$ I ,/y2+j(m-ri)$ H /v / 2 



h LtR = (m±n) 2 /4 



(12a) 
(12b) 



where we have fixed the compactification radius at r = 
V2. The labels to and n are related to charge and mo- 
mentum, respectively. The relevant operators that are 
constructed as products of these are listed in Table I. 
Note that the model is tuned to the Kosterlitz-Thouless 
(KT) transition, where the cosine term, cos(2$/-\/2), is 
precisely marginal (Hl + hn = 2). 

To see that this model is supersymmetric, remember 
that the left (right) moving supercharges have scaling 
dimension (h L ,h R ) = (3/2,0) {(h L ,h R ) = (0,3/2)). It 
follows that any c = 3/2 theory has an Af = (1,1) su- 
persymmetry, generated by the left- and right-moving su- 
percharges tpLd&L and iPrB^r. For the compactification 
radius r — V2 there is an additional Af = (2, 2) super- 
symmetry, generated by the two left- and right-moving 
supercharges ipL,Rexp[±iV2®Rx]- 

We now return to the lattice model and briefly review 
the argument that led to propose the c = 3/2 super- 
conformal field theory as its continuum theory and some 
results that are in agreement with this conjecture. We 
first note that the Witten index, W = k + 1, is equal 
to 3 for the 2nd minimal model. This is in agreement 
with the Witten index of the lattice model, given by 
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W = Tr(— 1)5^ which is 3 if the length is a multiple 
of 4. Furthermore, a typical configuration in the ground 
state is IL^l j)|0)> which we schematically depict 
as . . . OlOiOIO^OIofo . . . . This suggests the possibility 
of two types of domain walls associated with defects of 
charge ±1/2 (OHO), on the one hand, and with neutral 
defects (JOJ) on the other hand. As was proposed in 
Ref. 1 , it is suggestive to associate the Z 2 lattice parity 
symmetry that maps I into I and vice versa with the 
Ising sector in the c = 3/2 theory. The U(l) symmetry 
associated to the charge degrees of freedom then simply 
corresponds to the bosonic sector. Further evidence for 
this picture was given in Ref. 13, where two exactly sol- 
uble limits of this model were considered by perturbing 
away from the supersymmetric point. It was shown that 
in these limits, the model exhibits an Ising transition and 
a KT transition, respectively. It was suggested that these 
two transitions coincide at the supersymmetric point. We 
briefly review the soluble limits in App. A. Finally, it was 
shown in Ref. 4 that the supersymmetric model is gap- 
less by numerical studies of a boundary twist using exact 
diagonalization. These numerical studies, however, also 
revealed the strong finite-size effects that made it thus 
far impossible to conclusively identify the c — 3/2 theory 
as the continuum theory of this model. 



A. Marginal operators 

The operator content of our candidate theory in prin- 
ciple allows for a variety of marginal operators with scal- 
ing dimension h = Kl + Hr = 2. These could serve 
as an explanation of the strong finite-size effects ob- 
served in numerical calculations. On the other hand, one 
might naively expect that these operators are excluded 
by the explicit supersymmetry on the lattice. We have 
found, however, that there is a special combination of 
two marginal operators that preserves Af = 2 supersym- 
metry. The two marginal operators are the usual cosine 
term for the boson, cos(\/2<I>), and the more unusual op- 
erator (d x $)iipRipL- In Section VIII, we calculate the RG 
equations for the couplings to these operators and argue 
that they are not excluded by symmetries for a certain 
choice of the couplings. 

As a result of this, slowly decaying corrections to the 
energy and effective central charge are obtained (see Sec- 
tion VIII for more details). Corrections to the free- 
energy density due to marginal operators were calculated 
in Ref. 17. Given our candidate CFT, we find for the 
scaling of energy levels 



the central charge due to the presence of a marginally 
irrelevant operator is ' 



E(L) = E/L + ai /{Ly/log(L)) + 



(13) 



where ai is some non-universal parameter. Clearly, one 
has to go to extremely large system sizes to reliably ex- 
tract the universal value E = %l + Hr — c/12. For the 
central charge, however, it turns out that the finite size 
corrections are strongly suppressed. The scaling form of 



c(L) 



(14) 



where a is a fit parameter. We will discuss the impli- 
cations of this effective central charge for our numerical 
approach in the next section. 



IV. PHASE DIAGRAM 

Our candidate c = 3/2 CFT can be characterized as 
the product of a free boson (Luttinger liquid with c = 1) 
and a critical Ising model (c = 1/2). Furthermore, the 
bosonic part of the theory is tuned to the Kosterlitz- 
Thouless (KT) transition where the staggering operator 
is exactly marginal. Such a critical point is expected to 
appear as a multicritical point in a two-parameter phase 
diagram. One would expect four phases adjacent to this 
critical point, corresponding to ordered and disordered 
phases of the charge and the Ising degree of freedom. 

Our system has charge conservation, providing an ob- 
vious identification of the charge degree of freedom. Since 
we consider spinless fermions, the Ising part of the theory 
cannot be related to particle spin. As was proposed in 
Sec. Ill, we will see instead that it is related to the lat- 
tice parity symmetry of exchanging the upper and lower 
chains. We find that the phase diagram can most easily 
be obtained by tuning the following three parameters: i) 
the rung hopping t±, which favors states odd under ex- 
change of the sites on a rung, ii) a two-body density re- 
pulsion J (defined below), which favors charge ordering, 
and iii) the filling. While in general only two parameters 
should need to be tuned to obtain such a phase diagram, 
we have found that one of the four possible phases ad- 
jacent to the supersymmetric point cannot be realized 
without changing the filling. 

The perturbations can be written using the operators 
introduced in Section II in the following form: 



H, 



pert 



J Yl i n i,l n i+hl + n i,l n i+l,l) ( 15 ) 



Note that in these units, t± = J = corresponds to the 
supersymmetric point. At tj_ = —1, there is no rung 
hopping because the perturbing term exactly cancels the 
one in the original Hamiltonian. The term proportional 
to J corresponds to the potential term H V 2 of Eqn. (5); 
thus, at J = —2, this term is removed from the Hamil- 
tonian entirely. Unless otherwise mentioned, we will as- 
sume that the chemical potential has been adjusted such 
that the ground state remains at quarter filling. The fill- 
ing is enforced exactly by using quantum numbers in the 
DMRG calculations. In some phases, we tune the filling 
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FIG. 1. Phase diagram for the model in the t±-J plane. The 
labels show the shorthand notation for the phases: C corre- 
sponds to charge order, I to Ising order, and CI to charge and 
Ising order. The phases are discussed in more detail in Sec- 
tion IV. Our numerical approaches to establishing the phase 
boundaries are discussed in Section VI B. 




FIG. 2. Illustration of the two gapped phases. Top panel: 
charge ordered, Ising disordered phase. The ellipses indicate 
sites occupied by a fermion in the antisymmetric state (di+ — 
di-)\0). Bottom panel: charge and Ising ordered phase. The 
bigger red sites indicate occupied sites, the other sites are 
empty. 



away from the ground state filling N = L/2 by inserting 
a small number of holes. 

Perturbations of this model away from the supersym- 
metric point were also considered in Ref. 13. It was found 
that in the limit t± = — 1 the model can be solved exactly 
by mapping it to an XXZ Heisenberg chain. The KT 
transition between the charge ordered and charge disor- 
dered phase is found to be at J = in this limit. Further- 
more, in the limit J — > oo the model maps to an effective 
Ising model, with an Ising transition at t± = —1/2. We 
give more details on these soluble limits in App. A. 

We expect four possible phases adjacent to the su- 
persymmetric point, which are characterized by charge 
(dis-)order and Ising (dis-)order. We find that for the 
case without doping away from quarter filling, only three 
phases are realized: (i) an Ising ordered, charge ordered 
phase which is fully gapped, (ii) a Ising disordered, charge 
ordered phase, which is also fully gapped, and (iii) an 



Ising ordered, charge disordered Luttinger liquid phase. 
We will refer to these phases using the shorthand nota- 
tion C for charge ordered, I for Ising ordered, and CI for 
charge and Ising ordered. When the system is doped be- 
low quarter filling, the holes become itinerant in both of 
the fully gapped phases, giving rise to a gapless charge 
mode. This allows us to realize the (iv) Ising disordered, 
charge disordered phase. The phase diagram obtained 
without doping is shown in Fig. 1. Our methods for es- 
tablishing these phase boundaries will be discussed in 
Section VI B. 

We now discuss the exact ground states in two limits 
corresponding to the two gapped phases, providing an 
intuitive picture for these phases. In the C phase, we 
expect a state reminiscent of a valence-bond solid, where 
every other rung is occupied by an odd combination of 
the two states on the rung, i.e. the antibonding state. In 
the limit ij_ — >■ oo, the state should be a product state 

i^ = (i) L/2 (4 iS -< I )(< J -< I )...|o>. (i6) 

This state is shown in the upper panel of Fig. 2. In this 
limit, the translational symmetry is broken to a four-site 
unit cell and the state should be twofold degenerate. At 
finite values of t± , we expect this degeneracy to be lifted 
since the symmetric and antisymmetric combination of 
the two translated product states are separated by an 
exponentially small energy gap. 

The other fully gapped phase, CI, is characterized by 
particles alternating between the upper and lower chain 
with an empty rung between each particle. In the limit 
t± — > 0, J — > oo, the state is 

i^) = 4 )J 4 iI 4 )S ...|o). a?) 

This is shown in the lower panel of Fig. 2. Due to transla- 
tional symmetry breaking, one would expect an four-fold 
degeneracy in the extreme limit, which is again lifted at 
finite values of t± and J. 

The Ising ordered, charge disordered phase I can be 
understood by considering the CI phase discussed be- 
fore, and allowing an overall charge mode. In particular, 
due to Ising order, the density correlations still show the 
alternating structure between the upper and lower chain. 
Similarly, one would expect an Ising disordered, charge 
disordered phase upon doping the C phase, where the 
antibonding states of phase (i) become itinerant. How- 
ever, the infinite nearest-neighbor repulsion implies that 
two fermions cannot simultaneously be in the antibond- 
ing state on adjacent rungs. The Ising disordered, charge 
disordered phase can therefore not be realized using the 
parameters discussed so far, but instead must be realized 
by hole-doping the system away from quarter filling. 

Phase C and CI are separated by a line of c = 1/2 
Ising transitions, and phase I and CI are separated by a 
line of KT transitions. At t± = —1, the system can be 
solved exactly and the transition is found to be exactly 
at J = (see App. A). The line of KT transitions and 
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the line of Ising transitions join at the supersymmetric 
point. The nature of the transition from phase I to phase 
C is unclear and will be discussed below. 

When the system is doped away from quarter filling by 
inserting holes, these holes become itinerant in both of 
the fully gapped phases, giving rise to a gapless charge 
mode. The central charge in these phases therefore be- 
comes c = 1. The line of KT transition separating the 
two Ising ordered phases disappears in this case and only 
two phases remain. These are separated by a line of Ising 
transitions, which appear as a line of c = 3/2 transitions 
within a Luttinger liquid phase. 

V. METHODS 
A. DMRG 

Our calculations are carried out using the density ma- 
trix renormalization group method. 1D ~ 22 This method is 
generally formulated for Hilbert spaces with a tensor- 
product structure and therefore does not allow the im- 
plementation of the hard-core constraint on the level of 
the Hilbert space. Instead, we add a penalty term to the 
Hamiltonian that increases the energy of configurations 
with occupied nearest-neighbor sites. We find that this 
term does not have to be very strong to obtain reliable 
convergence as the hopping term is dressed with the pro- 
jection operator and therefore only acts on the allowed 
subspace. In addition, we choose a basis where the two 
sites on a rung are treated as a single site, allowing us to 
implement the hard-core constraint exactly for this rung. 

The approximation made in DMRG calculations can be 
systematically refined by increasing the number of states 
M kept in the renormalization procedure. The algorith- 
mic cost grows as C(Af 3 ), which limits the number of 
states to a few thousand in practical calculations. While 
for gapped systems, a bond dimension of a few hundred 
is generally sufficient independent of system size, criti- 
cal systems are more challenging because the number of 
states has to be increased as some polynomial of the sys- 
tem size when the thermodynamic limit is approached. 
In addition, periodic boundary conditions strongly in- 
crease the number of states needed when using a standard 
DMRG approach. While improved schemes exist, 21 ' 22 we 
use a well-tested and numerically robust standard ap- 
proach. While most of our simulations are carried out 
using M = 1000 states, we confirm results for up to 
M — 4800 states for some long periodic systems. This 
limits us to systems of length up to L = 100. 

B. Identification of the conformal field theory 

Correctly identifying the field theory that describes the 
continuum limit of a given microscopic model is a noto- 
riously difficult problem. In general, one has to resort 
to numerical simulations, which are usually restricted to 



finite systems. While the most commonly used method 
of identifying the CFT is by fitting the spectrum ob- 
tained with exact diagonalization for small systems to a 
spectrum obtained directly for the CFT, we use a differ- 
ent approach in this paper: we first establish the central 
charge at the supersymmetric point by studying the en- 
tanglement entropy of the system. We then move on to 
study the phase diagram of the model with various per- 
turbations, where we again rely on entanglement entropy 
- augmented with calculations of correlation functions or 
structure factors - to establish the phases. Finally, we 
study a perturbation where in one limit, the model in the 
exactly soluble case of a chain is recovered. 

As discussed above, finite-size corrections to the cen- 
tral charge due to marginal operators are strongly sup- 
pressed. Calculating the central charge therefore appears 
as a much more viable approach to identifying the CFT 
than studying the energy spectrum. We will now elabo- 
rate on how this can be done in numerical calculations, 
and then discuss possible corrections in finite systems. 

We consider the entanglement entropy S(p) = 
— Tr plogp where p is the reduced density matrix for 
some block of sites. If we denote the entanglement en- 
tropy between a block of I contiguous sites in an infinite 
system and the rest of the system as 5(0, we have for 
the entropy of this block 23 

5(1) ~ | log/. (18) 

Using this relation, the central charge can be extracted 
by calculating the entanglement entropy at the center of 
a finite system of length L and performing a fit to 

Tl C 

5(/ = L/2)~ylogL, (19) 

where n = 1 for open and n — 2 for periodic boundary 
conditions. For reasonably large systems, which can be 
simulated with the DMRG method, such a fit often gives 
accurate results for the central charge. 

For a finite block embedded in a finite systems, 
Eqn. (19) only holds approximately. A more appropri- 
ate relation for finite systems is obtained in Ref. 24. The 
entropy for a block of I sites in a finite periodic system 
of length L is 

5(/) = | log sin ^ +S , (20) 

and for I sites at the end of an open system we have 

S(Z) = £ log sin ^)+S . (21) 

In the constants Sq, we have summed up several universal 
and non-universal contributions, which are not relevant 
for our purposes. This expression opens up a second way 
of determining the central charge: instead of simulat- 
ing several different system sizes and performing a fit to 
Eqn. (19), one can simulate only a single system size and 
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FIG. 3. Finite-system gap for the ladder with open bound- FIG. 4. Effective central charge obtained for a ladder with 
ary conditions. Calculations were performed with M = 1000 antiperiodic boundary conditions using M = 4800 states, 
states. The dashed line shows a fit to Eqn. (23). 



calculate the entanglement entropy for various block sizes 
I. A fit to Eqns. (20) and (21) will then yield the desired 
central charge. We will show numerically that the cen- 
tral charge obtained in this way contains the corrections 
discussed in Sections III A and VIII. 

We note in passing that additional corrections occur 
for open boundary conditions. These were first observed 
numerically in Ref. 25 and later explained analytically 
in Refs. 26 and 27. The correction can be attributed to 
Friedel oscillations from the boundary, and is well fit by 



L nl 
— sin — 

7T L 



-K 



(22) 



where K is the Luttinger liquid parameter in the case 
of c = 1 CFTs, and related to the scaling dimension of 
relevant operators otherwise. The correction can heuris- 
tically be explained by a small dimerization on open lat- 
tices. Due to this correction, we will restrict our calcula- 
tions to periodic systems. 



VI. RESULTS 
A. The SUSY point 

We start by demonstrating that the system is indeed 
gapless with a dynamical critical exponent of z = 1, i.e. 
that A(L) ~ and since the correlation length is ex- 
pected to diverge at a critical point, A(L) ~ (to 
leading order). This is a the necessary condition for the 
system to be described by a conformal field theory in the 
continuum limit. 

In Fig. 3, the finite-size charge gap, defined as A(L) = 
E(N = L/2 + l)+E(N = L/2-1), where N is the num- 
ber of particles, is shown. Note that the ground state 



has iV = L/2, and E(N = L/2) = by supersymme- 
try. Since this calculation is performed for open bound- 
ary conditions, a modest number of states in the DMRG 
procedure is sufficient to obtain good accuracy even for 
systems of length L > 100. We perform an extrapolation 
by a fit to 



A(L) 



ax a 2 
L L v / logX 



(23) 



We expect in particular that A = 0, i.e. the sys- 
tem is gapless. The subleading terms are included due 
to the presence of marginal operators in the candidate 
CFT (see Sees. Ill A and VIII). Such corrections terms 
were discussed in Refs. 17 and 18. The good agreement 
with the expected scaling confirms that the system is in- 
deed described by a conformal field theory perturbed by 
marginal operators. The accuracy could be improved by 
including further higher-order corrections. 

We next move on to study the entanglement entropy 
at the supersymmetric point. In order to select a unique 
ground state and avoid oscillating terms from the bound- 
aries, we use antiperiodic boundary conditions for this 
calculation. Due to this choice of boundary conditions, 
we have to use a large number of states in a conven- 
tional DMRG procedure in order to obtain accurate re- 
sults also for large systems. In Fig. 4, data is shown that 
was obtained using M — 4800 states, which is sufficient 
to exhaust the entropy for these system sizes. The cen- 
tral charge we obtain shows large finite-size corrections, 
which is unusual. We will now show that these correc- 
tions can consistently be explained by the presence of 
marginal operators. 

To account for finite-size corrections due to marginal 
operators, we perform an extrapolation of the central 
charge c(L) measured in a finite system of length L us- 
ing Eqn. (14), where we use both c and a as fit parame- 
ters. We find good agreement and an extrapolated value 
reasonably close to the expected value c = 1.5. This 



gives strong indication that the central charge is indeed 
c = 3/2. It also shows that the corrections to the cen- 
tral charge originally obtained for the free energy density 
can analogously be applied to calculations based on the 
entanglement entropy. 



B. Phase diagram 

To establish the phase boundaries, we calculate the 
central charge for a fixed system size L = 24 on a fine 
mesh of points in the (t±,J) plane. The phase bound- 
aries are then easily extracted from the behavior of the 
central charge. Our results are shown in Fig. 1. As ex- 
pected, the I phase appears as an extended region with 
central charge c = 1. The precise location of the transi- 
tion from this phase into the CI phase is known to lie at 
J = for t± = and tj_ = -1. For -1 < t± < 0, the 
location is unknown, and due to the Kosterlitz-Thouless 
nature of the transition is difficult to determine numer- 
ically with a high accuracy. From our results for few 
system sizes, however, we can exclude significant devia- 
tions from the J = line. The transition from the CI 
to the C phase appears as a line of c = 1/2 transitions, 
as expected. This is expected to approach t± = 0.5 for 
J — > oo, which is consistent with our observations. Note 
that the numerical results for the line of Ising transitions 
agrees well with the large J result, t± = — l/2 + l/(J+2), 
even as J — > (see App. A). 

The nature of the transition from the I phase to the 
C phase is at this point unknown. Several scenarios are 
possible: i) A line of c = 3/2 transitions emanating from 
the supersymmetric point where the Ising transition and 
the charge ordering transition take place simultaneously. 
This requires fine-tuning and seems plausible only in the 
presence of an additional symmetry, ii) Without fine- 
tuning, the Ising transition could either be in the gapped 
or the charge-disordered phase. In the first case, it would 
correspond to a line of c = 1/2 transitions, whereas in 
the latter case c = 3/2 should be observed with c — 1 on 
either side. 

Numerically, we find no evidence of two separate tran- 
sitions. Instead, the central charge shows a single, yet 
broad peak. This persists to system sizes up to L ~ 70 
with antiperiodic boundary conditions. The position of 
the peak shifts significantly with system size. While this 
behavior seems more consistent with the first scenario, 
we observe an unusually large value of the central charge 
approach c = 2 for J — > —2. This cannot be explained 
in this scenario. Additionally, since the observed peak is 
very broad, one could argue in favor of the second sce- 
nario where the two transitions are so close to each other 
that they cannot be resolved on the accessible system 
sizes. In conclusion, the precise nature of this transition 
remains an open question. 

To further corroborate the nature of these phases, we 
study the density structure factor. Defining first the odd 




FIG. 5. Structure factor for L = 24 and antiperiodic bound- 
ary conditions. Top panel: Charge and Ising ordered phase, 
J = 1, t± = — 1. Second panel: Charge disordered, Ising 
ordered phase, J = — 1, t± = — 1. Third panel: Charge or- 
dered, Ising disordered phase, J — 0, t± = 0.65. Bottom 
panel: Supersymmetric point. 



and even density on each rung, 

we can define the even and odd structure factors 



(24) 



S±(k) = J2e iM («n? d ) ~ (n±)(n^ d )) (25) 

In a translationally invariant system, the index i can be 
fixed to 0. 

In Fig. 5, the numerically calculated structure factor 
for a fixed system size and four different values of the 
parameters, corresponding to the three different phases 
and the supersymmetric point, are shown. By inspection, 
Ising order corresponds to peaks in the S- structure fac- 
tor at k = ir/2 and k = 3ir/2, whereas charge order cor- 
responds to a peak in S+ at k = tt. To understand this 
behavior, one can consider the ground state in the limit- 
ing cases given in Eqns. (16), (17). In the C phase, n_ 
vanishes on all rungs, whereas n+ is large on every other 
rung. This leads to the peak at k = ir/2 and no other 
strong features, as shown in the third panel of Fig. 5. In 
the CI phase, particles sit on every other rung, so that 
n + is finite only for those rungs; n_ is also finite only 
on these rungs, but in addition oscillates in sign. There- 
fore, we find a peak in S+ at k = n and S- at k = n/2, 
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FIG. 6. Strength of the density correlation function (norij±) (shown only with j > 0) for values of the parameters representative 
of the three different phases at quarter filling. 



k = 3tt/2. As discussed above, the I phase is most eas- 
ily understood by introducing an overall charge degree 
of freedom in the CI phase, but keeping the alternating 
structure. Consequently, the peak in S+ is strongly re- 
duced compared to those in S-, cf. the second panel of 
Fig. 5. The bottom panel, corresponding to the super- 
symmetric point, illustrates the proximity of this point 
to the ordered phases. 

In Figure 6, the correlation functions (^o J^i j) 
and (n j n i j ) for antiperiodic boundary conditions are 
shown, where the diameter of the circles indicates the 
strength of the respective correlation. This corresponds 
to the density structure that would be obtained if transla- 
tional symmetry was broken. The upper panel shows the 
system in the CI phase. Clearly, the density structure is 
consistent with the expectations and reminiscent of what 
is shown in the lower panel of Fig. 2. The center panel 
shows the density correlation in the I phase. Clearly, 
density correlations are weaker in this phase, which is 
expected due to the gapless mode. The Ising order is nev- 
ertheless still discernible. Finally, the lower panel shows 
the system in the C phase. For sites sufficiently far away 
from the O'th rung, the symmetry between the upper and 
lower chain is restored, which is expected for the Ising 
disordered phase. Charge ordering is also clearly visible, 
with every other rung showing much higher occupation. 

Finally, we study the behavior of the system when 
doped away from quarter filling. In Fig. 7, we show a 
numerically obtained phase diagram for L = 24 with 0, 2 
and 4 holes. The phase diagram was established by mea- 
suring the central charge for a fixed system sizes L = 24 
on a fine parameter grid. The phase boundary can then 
be observed as a line with c = 1.5, whereas we have 
c = 1 in the rest of the phase diagram as all phases have 
at least a gapless charge mode. While the location of the 
phase boundaries is not of much interest here, this con- 
firms that the Ising and charge disordered phase can be 
realized in a wide parameter range by doping the system 
away from quarter filling. 

VII. PERTURBATION TO THE CHAIN 

The second supersymmetric minimal model contains 
a supersymmetry-preserving relevant operator that in- 
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FIG. 7. Location of the Ising transition in the t±-J phase 
diagram for doping away from ground state filling. The tran- 
sition is determined from the peak in the central charge as t± 
is varied for fixed J. Calculations are performed for L — 24 
with the number of particles N = L/2 — Nhoie- 
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FIG. 8. Central charge for systems perturbed away from the 
isotropic supersymmetric point by staggering on two sublat- 
tices, cf. Eqn. (27). The chemical potential [s, is only added 
on sites on sublattice Sa. 
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duces a flow towards the first supersymmetric minimal 
model with c = 1. We can therefore obtain further evi- 
dence for the conjectured continuum theory by identify- 
ing a microscopic perturbation that corresponds to this 
operator in the field theory. Such a perturbation should 
meet three criteria: i) In the limiting case, the ground 
state of the perturbed ladder system coincides with that 
of the chain, and in particular is critical with c = 1. ii) 
The Ising sector becomes gapped when turning on the 
perturbation, but the charge sector remains gapless and 
the unperturbed limit and the chain limit are connected 
by a line of critical systems with c = 1. iii) The operator 
should preserve supersymmetry on the lattice. 

Supersymmetry-preserving operators can be con- 
structed using the staggering discussed briefly in Sec- 
tion II. We find that for a specific choice of staggering, 
the limiting case indeed corresponds to the chain and 
the perturbation thus meets our first criterion. We will 
demonstrate numerically that the central charge remains 
at c = 1 between the chain and the ladder limits, which 
confirms the second criterion. 

The specific staggering divides the lattice into two sub- 
lattices Si and S2 with 

S 2 = {(An -3,1), (An - 1, I )|n = 1, . . . , L/4}, (26) 

where (i, I) ((i, I)) indicates the upper (lower) site of the 
i'th rung. The supercharge is then defined as 

ieSi i£S 2 

Note that S2 coincides with the occupied sites in the 
bottom panel of Fig. 2. 

For z = 1, the system coincides with the homogeneous 
ladder system studied before. In the limit z = 0, there 
is no hopping term between the two sublattices and the 
particle number on each sublattice becomes a good quan- 
tum number. One can easily determine that one possi- 
ble ground state of the system with periodic boundary 
conditions has all particles on sublattice S2, which is a 
non-entangled product state. Since the total three-fold 
ground state degeneracy must remain unchanged under 
staggering, there must be two ground states where all 
particles are on sublattice Si. This sector is equivalent 
to a chain with length L' — 3L/2, where L is the length 
of the ladder, and filling 1/3. The two-fold ground state 
degeneracy is then consistent with the degeneracy of the 
chain model, and we expect this state to be critical with 
c = 1 and described by the first TV = 2 supersymmet- 
ric minimal model. We can easily select the ground state 
sector corresponding to a critical chain by adding a chem- 
ical potential on sublattice S2 such that particles on that 
sublattice are penalized. For z > 0, the number of par- 
ticles on each sublattice is not a good quantum number 
and we expect the ground states to mix. 

We numerically calculate the central charge for values 
of z between the limiting cases z = and z = 1. Our 
results are shown in Fig. 8. The solid lines show results 



without chemical potential. Clearly, the central charge 
drops to c = 1 away from the z — 1 point. For small 
values of z, it becomes increasingly difficult to obtain the 
critical ground state because the DMRG calculation is bi- 
ased towards the low-entanglement solution of particles 
sitting only on the S2 lattice. While this is a true ground 
state only at z = 0, it is very competitive in energy also 
for z > 0. This is remedied by adding a chemical poten- 
tial on 5*2 for small values of z (dashed lines). This gives 
c = 1 also for small values of z. Our data are therefore 
consistent with a c = 1 phase for < z < 1, which is 
connected to the supersymmetric model on the chain. 

VIII. EFFECTIVE FIELD THEORY 

In this section we consider the continuum theory un- 
derlying the supersymmetric multicritical point and the 
role of marginal perturbations. The second superconfor- 
mal minimal model with c = 3/2 describes a compactified 
boson, $ = $1, + and a free Majorana fermion, t/>, 
with J\f = (3,3) supersymmetry. In principle the opera- 
tor content of this theory allows for a variety of marginal 
operators with scaling dimension h = hi, + fiR = 2, that 
could cause the strong finite size effects that we observe 
numerically. Most of these operators, however, are ex- 
cluded either by lattice symmetries or by supersymme- 
try. As mentioned in Section III A, we have found that 
there is a special combination of two marginal operators 
that preserves Af = 2 supersymmetry. The two marginal 
operators are the usual cosine term for the boson and the 
more unusual operator (d x ^)iipR^pL- In the following we 
analyze the renormalization group flows of the c = 3/2 
theory in the presence of these operators without impos- 
ing supersymmetry. At the end of this section we will 
discuss the flow along the direction that preserves super- 
symmetry to see the finite-size effects resulting from the 
presence of these marginal operators. 

The renormalization group flow of the c = 3/2 the- 
ory in the presence of the operator (d x ^>)iipRipL was 
worked out in Ref. 29 by considering Wilsonian one-loop 
momentum-shell RG. In order to analyze the cosine term 
it will be more convenient to work in real-space and use 
the Operator Product Expansion formulation of RG. 
We consider the lagrangian density C = C& + Cjr + £ int , 
where 

£ *=2^(>^ 2 + ^ 2 ) 

= ipR(d T + -d x )ip R + i) L (dr 9 x )ipL 

1 1 

Ant = -A(a x $)*VjW>z +0cos[\/2$]. (28) 

Here, v,K are the velocity and Luttinger parameter of 
the bosonic field, u is the fermion velocity, and A, g are 
coupling constants. Note that with these conventions 
the scaling dimension of the cosine term is h = K/2, 
which is marginal at K = A. We treat the interaction 
terms perturbatively. In a real-space RG approach one 
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takes care of UV divergencies by introducing a real-space 
cutoff, a, which defines the minimal distance that two 
operators can approach each other. In an RG step this 
cutoff is increased, a — > (1 + 5£)a, while the effective 
action is kept fixed by renormalizing the couplings. We 
work in a scheme where we keep the bosonic velocity v, 
the normalization of the field $ and the unit prefactor 
of the term ^2 i=L R ipid T ipi fixed. To achieve this we 
allow the fermion field to rescale, ip — > ^/Zpip, and we 
introduce an anomalous dynamical exponent, z. Having 
z ^ 1 reflects the fact that the A-term breaks Lorentz 
invariance. For the field rescaling we find 



Z, 



f 



X 2 K 



8(v 



and for the dynamical exponent we obtain 



z = 1 



X 2 K 
I6uv 



The fermionic velocity, u, the Luttinger parameter, K, 
and the interactions A and g flow under RG. To second 
order in A and g we obtain the following RG equations 



du 
~dl. 
d\ 

M 
dK 
111 

dg 



iX 2 K 



1 



.(v + u) 2 



Auv / 



= 



K 



( 9 2 n 
V 2v 2 
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The first two equations and the dynamical exponent pre- 
cisely agree with the results of Ref. 29 (up to slight differ- 
ences in conventions 51 ). The latter two equations reduce 
to the Kosterlitz equations, when we set A to zero. In- 
deed, for A = the boson and Ising sectors are completely 
decoupled. Finally, we find that for u = v and K = 4 
there is a line given by 

A = 2V2vr 5 , 

where the RG equations are all zero to second order in the 
couplings. Note that for this value of K the cosine term is 
indeed marginal. It can be shown that precisely this line 
preserves Af — 2 supersymmetry. To prove this is beyond 
the scope of this paper and will be deferred to Ref. 32. It 
follows that coupling to this special combination of these 
two marginal operators is not excluded by the explicit 
lattice supersymmetry. 

The line of fixed points is likely to disappear at higher 
order in perturbation theory 2 ' ly '' 2 Here we argue that to 
higher order there will be a flow towards the fixed point 
A = g = based on the following general considerations. 
The central charge of the Af = 2 superconformal minimal 
models is given by c = 3k/ (k + 2) = 1, 3/2, 9/5, 2, . . . . 
Our numerical analysis of the central charge gives a value 
close enough to c = 3/2 to be able to exclude the other 



minimal models. Finally, using the fact that for L — An 
the lattice model has 3 zero energy states in the Ramond 
sector (pbc) and 1 negative energy state in the NS sector 
(apbc), we can exclude other c — 3/2 theories with higher 
supersymmetry. 12 

As already mentioned in Section III A, the extremely 
slow flow along the supersymmetric line parametrized 
by A leads to strong finite-size corrections. Since there 
is no flow to second order, we expect that dXjdl = 
BX 3 + 0(A 4 ). The explicit computation of B will be 
presented elsewhere.' 2 For now, we will assume that 
B < 0, consistently with the DMRG results that in- 
dicate that the c = 3/2 fixed point is stable in the 
presence of supersymmetry. Solving the flow equation 
for A, we then find X(£) ~ 1/VI- For the scaling of 
energy levels this leads to a finite s ize scaling of the 
form, E{L) = E/L + at/(Ly/log(L)) + . . . , where at 
is a fit parameter. Clearly, one has to go to extremely 
large system sizes to reliably extract the universal value 
E = hi, + hn — c/12. For the central charge, however, 
it turns out that the finite size corrections are strongly 
suppressed. The scaling form of the central charge due 
to the presence of a marginally irrelevant operator is (cf. 
Eq. (14)) 17 . 18 



c(L) = c + a 2 A(L) 3 
1 



(yi0g(Z))3 



(29) 



where a is a fit parameter. 



IX. MICROSCOPIC DERIVATION OF THE 
FIELD THEORY 

In this section, we present a microscopic way to un- 
derstand the origin of the effective field theory (Eq. 28) . 
This provides a way to relate the physical observables 
to the continuum fields tp. In this section we define 
the axes such that x lies along the long direction of the 
ladder and y along the short direction. 

The microscopic Hamiltonian, Eq. 4, describes a one- 
dimensional system of particles with strong repulsive in- 
teractions. Due to the infinite nearest-neighbor repul- 
sion, there can be at most one particle per rung of the 
ladder; therefore, one can describe the system in terms of 
a "charge" (longitudinal) degree of freedom, correspond- 
ing to the position of the particles along the ladder, and 
an internal "spin" variable denoting the leg index of each 
particle. The symmetry of the system under reflection, 
y — > —y, implies that the spin degree of freedom has 
Ising-like (Z2) symmetry. 

For a quarter-filled chain, the ground state of the po- 
tential energy terms, Eq. (5)-(7), is a state in which ev- 
ery other rung is occupied by a particle whose position 
is staggered in the transverse direction to the ladder (as 
in the bottom panel of Fig. 2). Therefore, the internal 
Ising-like degree of freedom tends to order antiferromag- 
netically. Adding quantum fluctuations has two effects: 
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it introduces kink-anti kink pairs into the transverse Ising 
ordering pattern, and causes deviations from perfect crys- 
talline order along the ladder direction. Introducing a 
continuum Ising variable <r{x) and a displacement field 
$(x), one can write the density of the fluctuating crystal 
as 



n(x,y) = 



1 



2-s/27r 
A2 cos 



d x $ + Ai<t(x) cos 



2a/2 



Q 



2Q 



(30) 



where r = (a;, y) (with integer x, y) is the position vector 
along the ladder, Q = 27r(l/4, 1/2) is the wavevector of 
the perfect crystal, A12 are non-universal amplitudes, 
and the dots denote higher harmonics. Note that this 
periodicity is observed in the structure factors measured 
in the DMRG calculations (see Fig. 5). 

In the low-energy limit, we can deduce the most rel- 
evant terms in the Hamiltonian of $, a from symme- 
try considerations. The standard arguments for a one- 
dimensional fluctuating crystal lead to a free boson 
Hamiltonian for 4>. The lowest-order harmonic of the i> 
field which carries zero momentum (modulo a reciprocal 
lattice vector) is cos(v / 2 < I ) ); therefore this term can ap- 
pear in the Hamiltonian. The Z2 symmetry prohibits any 
odd power of the Ising order parameter a from appear- 
ing in the Hamiltonian. If the Ising field is tuned close 
to its critical point, it can be represented in terms of a 
free Majorana field ip, where the mass of the Majorana 
fermions iipRipL (or, equivalently, the energy operator of 
the Ising model) is tuned to zero. Finally, the product 
of the fermion mass and the boson density, lip 'Rip l9^> x , 
is allowed by symmetry, and can appear in the Hamilto- 
nian. 

The identification of the low-energy modes as long- 
wavelength fluctuations of charge density wave and Ising 
order parameters is supported by the proximity of charge 
and Ising ordered phases to the supersymmetric point, 
which are stabilized when supersymmetry is explicitly 
broken, as discussed in Sec. IV (see Fig. 1). 



X. CONCLUSION 

We have studied in great detail a critical supersym- 
metric lattice model. Despite being amenable to many 
analytic approaches, the continuum theory had not been 
firmly established prior to our work. This is mostly due 
to strong finite-size effects, which are caused by an un- 
usual marginal operator present in the superconformal 
field theory. The marginal operator breaks Lorentz in- 
variance, but preserves part of the supersymmetry. We 
demonstrate that the finite-size corrections are strongly 
suppressed in the entanglement entropy. This allows us 
to use entanglement entropy as a powerful probe for the 
system's properties. In particular, we show with care- 
ful finite-size extrapolation that the central charge of the 



model is c = 1.5, which is the value expected for the 
second superconformal minimal model. 

Augmenting entanglement entropy results with calcu- 
lations of correlation functions, we are able to establish 
a three-parameter phase diagram realizing four differ- 
ent phases adjacent to the supersymmetric point. The 
supersymmetric point appears as multicritical point in 
this phase diagram. We can relate the adjacent phases 
to the Ising and charge sector present in the candi- 
date superconformal field theory, providing additional 
strongly supportive evidence. Finally, we study a par- 
ticular supersymmctry-prcscrving perturbation which es- 
tablishes a relation of our model to the well-understood 
case of the chain, which is described by the first super- 
conformal minimal model. 

Our results underline that the entanglement entropy is 
a powerful probe for many-body systems since it can eas- 
ily be obtained using the density-matrix renormalization 
group and since it shows very robust features. 

The c=3/2 superconformal field theory, which we iden- 
tify as the continuum theory, is also found as the contin- 
uum theory of related supersymmetric models."' One of 
these models 11 exhibits both the J\f = 1 and the Af = 2 
supersymmetry of the continuum theory on the lattice. 
Therefore coupling to the marginal operator that we have 
found is excluded in this model. It would be interesting to 
see if relating this model to our model provides a way of 
identifying the marginal operator on the lattice. Ideally, 
such an identification would allow one to tune away the 
finite-size effects due to the marginal operator. 1 More 
generally, it might be possible to obtain qualitative veri- 
fications of the RG equations that we derived by study- 
ing supersymmetry preserving perturbations of the lat- 
tice model. 

Finally, we note that it has been suggested that mul- 
ticriticality is a generic feature of the supersymmetric 
lattice models and that the extensive ground state en- 
tropy in two dimensional systems is indicative of this 
feature. ■ The results of our paper may be seen as a 
small step in the very ambitious program of identifying 
multicriticality in these two dimensional models. An ob- 
vious next step is to carry out a similar analysis and 
explore the phase diagram by perturbing away from the 
supersymmetric point for a one dimensional model with 
an extensive ground state entropy such as the zig-zag 
ladder. 4 ' 12 
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Appendix A: Exactly soluble limits 

Here we review the results of the unpublished work 
of Ref. on the two soluble limits of the lattice model 
for hardcore spinless fermions on a 2-leg ladder. The 
starting point is the supersymmetric model on the square 
ladder given by the Hamiltonian in Eqn. (4). To this 
Hamiltonian the perturbations t± and J given in Eqn. 
(15) are added. The soluble limits are J — > oo, t± > — 1 
and t± = —1, J > —2, respectively. In the following we 
assume that we tune the chemical potential such that the 
ground state is always at quarter filling. 

Let us first consider the limit J — > oo , where we recover 
an Ising model. Large J implies a large energy penalty for 
particles to be on adjacent rungs. At quarter filling we 
can thus consider the low-energy subspace spanned by 
the configurations with a particle on every other rung. 
This subspace falls apart into two disconnected sectors 
with particles on all even/odd rungs. Projecting out the 
empty rungs, we are left with a Ising degree of freedom 
at each site representing a particle on the upper (lower) 
leg as f (I)- The hopping along the rung directly trans- 
lates into a spin flip. The 3-body term is always trivial in 
this low energy subspace, whereas the 2-body repulsion 
for two particles to be two sites apart on the same leg 
simply translates into a ferromagnetic repulsion between 
two neighboring Ising spins. Finally, there is another 
contribution to this ferromagnetic interaction from sec- 
ond order perturbation theory, where a particle virtually 
hops to a neighboring rung and back. Putting all this 
together the effective Hamiltonian for large J reads 

i/2 r 11 i 

i— 1 



where a x and a z are the usual Pauli matrices. It follows 
that for large J there is an Ising transition at t± = — ^ + 

2T7--V2. 

We now turn to the limit tj_ = — 1. In this limit there is 
no rung hopping, and since two particles cannot hop past 
each other the Hilbert space falls apart into many discon- 
nected sectors. From the large J limit, however, we know 
that for t± smaller than the critical value, the system is 
in an Ising ordered phase. In the Ising ordered phase the 
ground state is dominated by the configuration in which 
a particle on the upper leg is followed by a particle on 
the lower leg and vice versa. It follows that for t± = — 1 
the ground state must be in the sector where particles 
on the upper and lower leg alternate. In this sector the 
attractive 3-body term and the repulsive 2-body term 
for particles on the same leg add up to zero for all con- 
figurations. The only remaining interaction term is the 
2-body repulsion diagonally across a plaquette between 
a particle on the upper and a particle on the lower leg. 
It follows that for the effective Hamiltonian, we do not 
have to distinguish between particles on the upper and 
lower leg. Using a Jordan- Wigner transformation we can 
map the configurations of occupied and empty rungs to 
spin configurations. The nearest-neighbor hopping along 
the legs then translates into a nearest-neighbor spin ex- 
change term and the 2-body repulsion into a spin-spin 
interaction. The effective Hamiltonian for tj_ = — 1 reads 



r 2 + J 

H tx =-i — 4^ |^S' 4 + 5, i+1 + S. t S+ +1 H - SfS* +1 



where S z = a z /2 and S ± = {<r x ± i<j v )/A. The effective 
Hamiltonian is the well-known XXZ Heisenberg model. 
The continuum theory is the free boson. The KT tran- 
sition between the charge ordered and charge disordered 
phase coincides with the SU(2) symmetric point. Since 
this symmetry is also present in the lattice model the 
critical point is easily identified to be at J = 0. In prin- 
ciple one can consider perturbing away form the t± = — 1 
limit, but we will not do so here. 1 
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